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Abstract 

We investigate the electron population dynamics in an ensemble of 
nearly isolated insulating nanoparticles, each nanoparticle modeled as an 
electronic two- level system coupled to a single vibrational mode. We find 
that at short times the ensemble-averaged excited-state population oscil- 
lates but has a decaying envelope. At long times, the oscillations become 
purely sinusoidal about a "plateau" population, with a frequency deter- 
mined by the electron-phonon interaction strength, and with an envelope 
that decays algebraically as t~^^^ . We use this theory to predict electron- 
phonon dynamics in an ensemble of Y2 03:Eu^^ nanoparticles. 

1 Introduction 

Nanocrystals exist in a size regime which lies between that of atomic and bulk 
matter, thus making them ideal for the study of extreme quantum effects in 
condensed matter systems. In particular, the vibrational properties of nanopar- 
ticles are strikingly different than their bulk counterparts: A spherical nanopar- 
ticle of diameter d cannot support internal vibrations at frequencies less than 
about 2TTv/d, where d is a characteristic bulk sound velocity. Any property of 
the nanoparticle that depends on the vibrational spectrum, such as its thermo- 
dynamic properties or electron-phonon dynamics, will be very different at low 
energies than in bulk crystals. This will be especially true for nanoparticles — for 
example, in powder form — only weakly coupled to their surroundings. 

One way to probe the vibrational spectrum of a nanoparticle is to optically 
excite an electron-hole pair and study the intraband electronic energy relax- 
ation prior to radiative recombination j^]. However, the excitonic states, 
being only weakly localized, will suffer significant quantum-confinement effects 
in the nanoparticle, making comparison with bulk relaxation rather indirect. 
An alternative probe of the vibrational spectrum is provided by well-localized 
electronic impurity states in a doped nanocrystal p, H . The impurity states can 
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be used to probe both energy relaxation by phonon emission ||] and phonon- 
induced dephasing [Q. In these cases, the difference between the nanocrystal 
and bulk behavior is almost entirely a consequence of their differing vibrational 
modes. 

In a recent experiment, Yang et al. ^ measured the phonon emission rate 
in Eu-doped Y2O3 nanoparticles between two electronic states separated by 
3 cm~^ in energy. The experiment was performed on a powder of nanoparticles, 
prepared by gas-phase condensation, with a mean diameter of 13 nm. Regarding 
the nanoparticles as isotropic elastic spheres [H, the lowest internal vibrational 
mode (a five-fold degenerate torsional mode) has a frequency of about 11 cm~^. 
The excited-state population was found to decay exponentially with a lifetime 
of 27 fj,s, more than two orders of magnitude longer than that between the same 
levels in the bulk (about 220 ns). Although there is no quantitative theory 
available yet to explain their results, the inhibited phonon emission is consis- 
tent with a large suppression of the low-frequency vibrational density-of-states 
(DOS) expected for these small crystals. 

Experiments such as the one by Yang et al. ^ present an exciting opportu- 
nity to study the crossover of phonon dynamics from bulk to nanoscale systems. 
However, the following question immediately arises: Given the large apparent 
modification of the DOS by finite-size effects, is it still correct to expect expo- 
nential relaxation and to use perturbation theory (Fermi's golden rule and so 
on) to relate the relaxation rate to the phonon DOS? After all, the energy in a 
nanoparticle that is completely isolated from its surroundings would have to be 
exchanged between the electron and vibrational degrees of freedom in a Rabi-like 
manner, and no relaxation would be observed. Indeed, an isolated nanoparticle 
may be regarded as a phonon analog of a two-level atom in a cavity, which is 
known to exhibit oscillatory population dynamics Q . 

It should be possible to detect oscillatory population dynamics experimen- 
tally using the techniques of Ref . , but what would be the effect of the un- 
avoidable distribution of nanoparticle sizes and corresponding vibrational-mode 
frequencies? In the limit where each nanoparticle has a two-level system inter- 
acting with only the lowest-frequency vibrational mode, called the Lamb mode, 
each nanoparticle would exhibit vacuum Rabi oscillations |Q. However, the 
Rabi frequency, which is a function of the electron-phonon interaction strength 
and the energy mismatch between the two-level system and Lamb mode, would 
vary from nanoparticle to nanoparticle. 

The purpose of this paper is to investigate the electron-phonon dynamics of 
an ensemble of nearly isolated nanoparticles. We find that at short times the 
ensemble-averaged excited-state population oscillates in a fashion that depends 
on the detailed size distribution of the nanoparticles, but that at long times the 
oscillations become purely sinusoidal about a "plateau" population fraction, 
with the frequency of the oscillation determined by the electron-phonon inter- 
action strength alone, and with an envelope that decays algebraically as 
In the infinite-time limit, the excited-state population approaches the constant, 
finite plateau value. (Of course, the small damping of the Lamb mode, produced 
by the weak but nonzero interaction of a nanoparticle with its environment, will 
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eventually cause the electrons to relax irreversibly.) 

In the next section we derive a general expression for the average excited- 
state population in an ensemble of nanoparticles at zero temperature, and study 
its short-time and long-time behavior. The rotating-wave-approximation (RWA) 
we use there limits the application of our results to situations where the detuning 
energies (the energy mismatch between the two-level system and Lamb mode) 
are not much larger than the electron-phonon interaction strength. In Section 3 
we apply this theory to an ensemble of Y203:Eu'^"'" nanoparticles similar to that 
of Yang et al. ^ , and give quantitative predictions for the average excited-state 
population as a function of the mean nanoparticle size and standard deviation. 
Section 4 contains our conclusions and a discussion of questions for future in- 
vestigation. 

2 Ensemble-averaged population dynamics 

To begin, we consider a single doped nanoparticle. We assume this nanoparticle 
to be an isotropic elastic sphere with a single localized electronic two-level sys- 
tem embedded within. We limit our investigation to small detuning energies, 
thus allowing for the use of the RWA, which discards non-energy-conserving 
terms in the Hamiltonian |^. In addition, we will neglect the five- fold degener- 
acy of the lowest vibrational mode that would be present in a perfectly spherical 
nanoparticle and will assume a single nondegenerate Lamb mode. We shall 
return to this point below in Section 4. 

"Dissipation terms" in the Hamiltonian are neglected as well. By dissipation 
terms we mean interactions of the nanoparticle with its surroundings that allow 
vibrational energy to be carried away irreversibly. For weak dissipation, these 
terms would cause an exponential decay of the population at large times, causing 
the "plateau" population to fall exponentially to zero. This effect would be 
negligible at small times and therefore would have no impact on the short-time 
calculations given below. For large dissipation, the exponential decay would be 
noticeable at all times and would therefore render our results invalid. 

We assume that there is only one vibrational mode available for the electronic 
system to couple to — the next vibrational mode being so high in energy as 
to make the effect of its coupling negligible. Therefore, as stated, our doped 
nanoparticle system consists of a two-level atom coupled to single vibrational 
mode. The Hamiltonian, in units where ?i = 1, is given by 

H = ^ tac\ca -I- CJoa^a -I- ^ gaotlC^aCalip. + O"), (1) 

where a — 1,2. The first term in H is the Hamiltonian for a nonintcracting two- 
level system with energies t\ and £2, and fermionic creation and annihilation 
operators cj, and Cq. The second term is the Hamiltonian for a vibrational mode 
with frequency and phonon creation and annihilation operators and a. The 
third term is the ordinary first-order interaction between the two-level system 
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and the vibrational mode, gaai is the electron-phonon interaction energy, with 
5i2 = 521 = g, and the other terms equal to zero. 

The wave function can be written as a superposition of the two electronic 
states, together with all possible populations of the single vibrational mode, 

m)) = E C«„(i)e-^(^°+""«)* \an) , (2) 

an 

where 

1 («t)vt |0). (3) 
The coefhcients Can{t) satisfy the coupled equations. 



dtC2n{t) + zgy^e-'("°-^')*Ci,„+i(t) + zgV^e^(""+^^)*Ci,„_i(t) = 0, (5) 

where Ae = €2 — is the electronic energy-level separation. 

In the RWA, which is valid near resonance (defined by the condition ujq = 
Ae), these coupled differential equations reduce to 

dtCmit) + igV^e"''C2.n-i{t) = 0, (6) 

and 

dtC2n (t) + %g^W+\ e-'^'Ci.n+i (<) = 0, (7) 

where 

1/ = Wo - Ae (8) 

is the detuning frequency. 

These equations can be solved by Laplace transformation using the boundary 
condition CcmCO) = ^a2^na- The amplitude for the upper state is 

(9) 

where 



C2n{t) = SnO 



COS 



n = ^/^y^Tig^ (10) 

plays the role of the Rabi frequency in this problem 

The probability that the electron is in state a, irrespective of the number of 
phonons present, is 

No.it) =Y,\Can{t)\\ (11) 
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Then N2{t) = \C2o{t)\^ is given by 



and 



= 1 - ('^^) sin^ (yM > (12) 



Ni{t) = l-N2{t). (13) 



The dependence of Ni{t) and N2{t) on i/ will be suppressed for simplicity. 

In an ensemble of such nanoparticles, a variation in diameter yields a similar 
variation in detuning frequency, i'. We assume this to be a Gaussian distri- 
bution. To obtain values for the mean detuning j7, and standard deviation in 
detuning cr, we assume each nanoparticle to be an isotropic elastic continuum 
with stress- free boundary conditions. Then, as shown by Lamb the lowest 
vibrational frequency is 

Md) - (14) 

where Vt is the bulk transverse sound velocity and d is the nanoparticle diameter. 
The mean detuning is then 

T7 = wo(d)-Ae, (15) 

where d is the mean particle diameter. Using these relationships, we assume a 
distribution in detuning frequency given by 

P{v) ^ ^ . (16) 

The ensemble-averaged population of electronic state \a) is then given by 

/oo 
dvP{v)N^{t). (17) 
-oo 

The behavior of N2 (t) at short times is sinusoidal with a decaying envelope 
dependent on the specifics of the particle size distribution. We will give examples 
of the short-time dynamics below. 

The long-time behavior of N2 (t) can be obtained analytically by an asymp- 
totic expansion. We begin by writing ( [l^ ) as 

7V2(t)=iV2(oo) + ;^Re/(t), (18) 



where 



is a constant independent of time, and 



7V2(oo)^l--^ / d^-—-— (19) 
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Q2 , „2 e^'^^^\ (20) 

where 

ilrcs = 2.g (21) 

is the resonant Rabi frequency. It wiU turn out that the constant Ni(qo) is 
simply the value of f^iif) in the t — > oo limit. Values of this plateau population 
are given in Table 1. 
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Table 1: Plateau population values for as f — > oo. fires = 2g( is the 

resonant Rabi frequency. 



The integral l(t) can be evaluated at large times by analytically continuing 
V into the complex plane and expanding around the saddle-point at = 0. This 
leads to the asymptotic result 



7(t) « e-^'/"' A e'(*"-+"/4), (22) 



and hence 



+ «^oc,. (23) 

V 2 (T/iZrcs V"2rcs 

Note that at long times the population oscillates at the resonant Rabi fre- 
quency, independent of the mean detuning F. This occurs because the higher- 
frequency components tend to average out faster. However, the amplitude of 
the asymptotic oscillations decreases with 77. 

To illustrate the short-time behavior of N2{t) as a function of V and ct, we 
present several plots showing N2{t) as a function of tf^ros- In choosing values 
for V and cr, we attempted to cover the entire range of these parameters which 
could be reasonably addressed in the RWA. 

In Fig. 1, we study the resonant case where v = 0. By choosing three values 
of cr we are able to directly observe the effects that it has on the envelope 
function. As a increases, the envelope function decays faster and the time at 
which asymptotic behavior becomes observable decreases. In this case of F = 0, 
changing a has little effect on the frequency at which N2{t) oscillates. 

Fig. 2, where 77 = fires, illustrates the effects of an intermediate value of 77. 
The short-time envelope function decays faster and the plateau population is 
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Figure 1: Electronic population in state |2) for i' = 0.0. 



increased. As in Fig. 1, increasing a increases the rate at which the short-time 
envelope function decreases, but here it can be seen to lower the frequency of 
oscillation in N2{t). Unlike Fig. 1, Fig. 2 clearly shows how increasing a up to 
r^rcs increases the amplitude of oscillations at large times [see Eqn. (23)]. 
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Figure 2: Electronic population in state |2) for u = fl^ 
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The last plot, Fig. 3, illustrates the limits of the RWA. The effects of in- 
creasing v are now taken to the extreme: the short-time envelope functions 
decay very rapidly, making asymptotic behavior apparent at early times, and 
the plateau population approaches unity. Altering a has similar effects as in 
Figs. 1 and 2, only more pronounced. [A'^2(i) for cr/riros = 0.1 levels off to its 
plateau population shortly past iJli-os = 50.] 




Figure 3: Electronic population in state |2) for t^/rircs = 2.0. 



3 Population dynamics in Y203:Eu^+ nanoparti- 
cles 

Here we use the results of Section 2 to address possible future experiments. We 
examine the case of Y203:Eu'^+ nanoparticles, specifically the electron dynamics 
between two crystal-field split ^Di levels separated by Ae = 3 cm"-*^, such as 
in the experiment by Yang et al. To make quantitative predictions for 
such nanoparticles, the only parameters that need to be specified are the mean 
detuning frequency F, and the standard deviation in detuning a. We also need 
the electron-phonon interaction strength g. Our treatment requires that V be 
not much greater than fJi-es- 

In the case of Y203:Eu'^+ nanoparticles, we can estimate the electon-phonon 
interaction strength by using the experimentally observed bulk phonon emission 
rate |Q. We obtain a value for ^micron, the electron-phonon interaction strength 
in micron-size particles, and then scale this value to get g for the nanoparticles. 
The scaling is achieved by assuming that g varies with energy and system volume 
as it does in a bulk crystal. In principle, this method is only correct between 
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two systems with continuous vibrational spectra. However, using this technique 
should provide a reasonable estimate of g in the nanoparticle. 

In the micron-scale crystal, Fermi's golden rule states that for an electron- 
phonon interaction of the form given in Eqn. (|l|), the phonon emission rate 
is 

T^' = ^g'nAe)V, (24) 

where V is the crystal volume and Ae is the electronic level spacing. r(e), the 
phonon DOS (per volume) as a function of phonon energy e, is defined by 

n 

where cUn the phonon frequencies, and in a bulk crystal at low energies 
where 

/ \ -1/3 

is a branch-averaged sound velocity. 

The lifetime in the micron-sized crystal was observed to be t = 220 ns 
and in Y2O3 the sound velocities are approximately ||^, 

= 6.7 X 10^cms-\ (28) 
=4.3 X 10^cms^\ (29) 



and 



g is then given by 



4.7 X 10^ cm s"^ (30) 



The mean diameter of these nanoparticles was about 5 /im, leading to 

ffmicron = 2.04 X lO"*^ Cm'^. 

In a bulk crystal (with deformation-potential coupling), g varies as (e/y)^/^, 
so we scale ^micron in e and V to get a value for gnano, the electron-phonon 
interaction strength in the nanoparticle, 

-in — 1 \ 1/-^ /io \ —3/2 

10 cm \ / 13 nm ^ 



o -1 ; I r ; .9microni {^^) 

3 cm ^ y \ 5 fim J 
2.8 X 10-2 cm"\ (33) 
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Then the resonant Rabi frequency, as defined in Eqn. is about 

r2,es = 5.6 X 10^2 cm-i. (34) 

In the experiment of Yang et ai, the mean particle diameter was 13 nm 
and the standard deviation in particle size was 5 nm By Eqns. (|l^) and 
(p^, these values yield V = TlJlj-os and a = 4417i-es- This detuning is large 
enough that higher-energy vibrational modes can no longer be ignored, making 
our single mode treatment inapplicable. Furthermore, these values are too far 
off resonance to be reliably addressed by the RWA. 

To observe the behavior predicted in Section 2 for Y2 03:Eu'^+ nanoparticles, 
the mean diameter should be about 46 nm, and standard deviation in diameter 
no more than about 3 nm [see Eqns. (|l|) and 



4 Discussion 

We have shown that the unavoidable size-dispersion in a collection of nanopar- 
ticles can effect the ensemble-averaged electronic population dynamics consid- 
erably, affecting even the oscillation frequency at long times. However, our 
analysis has been restricted to near-resonant conditions in an idealized single- 
mode nanocrystal. A comprehensive theory of electron-phonon dynamics in 
nearly isolated nanoparticles will have to address a number of additional issues. 

First, the vibrational modes of a nanoparticle will be broadened due to 
intrinsic mechanisms and from interaction to the environment. For example, 
anharmonicity will broaden the Lamb mode at finite temperature as will 
mechanical coupling to a substrate or to a cluster of other nanoparticles pOf . 
This broadening will cause the envelope of the population oscillations to decay 
exponentially at long times. For weak damping, however, we do not expect any 
qualitative changes in the short-time dynamics. 

Second, nanoparticles will generally have large detuning frequencies, making 
the inclusion of higher- frequency modes necessary, and making the RWA invalid. 
One way to include these effects would be to do exact-diagonalization studies 
in models with a truncated Hilbert space; for example, having bounded phonon 
occupation numbers for one or more modes. 

Finally, it will be necessary to understand the effects of vibrational-mode de- 
generacy, which has been neglected here. The presence of degeneracy means that 
the Lamb mode has a vector-like internal degree of freedom, the components 
of the vector describing the phonon amplitude in each branch. Although there 
has been work done in the quantum optics literature on the multi-mode gener- 
alizations of the our Hamiltonian rtl ) . that work has focused on the two-photon 
resonance case where Ae = 2ujo ^1, |l^. One way to approach the degenerate 
case nonperturbatively would be to assume N degenerate modes in the limit 
of large N. By analogy with other quantum systems with internal degrees of 
freedom, like particles with N spin components, it is reasonable to expect that 
quantum effects will be diminished in the large N limit. In a nanoparticle this 
would suggest a suppression of the Rabi oscillations. 
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